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Abstract 

Renormalization group transformations for Schrodinger equation are performed in 
(/j 4 and in Yang- Mills theories. The dependence of the ground state wave functional 
on rapidly oscillating fields is found. For Yang-Mills theory, this dependence restricts 
a possible form of variational ansatze compatible with asymptotic freedom. 



1 Introduction 



It is widely accepted that many physical properties of gauge theories in the confining phase 
are determined by the structure of their vacua. In the literal sense of the term, the vacuum is 
described by a ground state wave functional - the lowest energy solution of the Schrodinger 
equation. To find a complete solution of the Schrodinger equation in Yang-Mills theory is 
by no means feasible, but the asymptotic freedom should allow to determine the dependence 
of the ground state wave functional on rapidly oscillating fields, since this dependence is 
responsible for the renormalization. 

General properties of the renormalization of functional Schrodinger equation were stud- 
ied within the perturbation theory. It was shown that both the Hamiltonian and the wave 
functional are renormalizable by usual counterterms up to multiplicative redefinition of the 
field variables [jlj. Practically, the most convenient way to perform renormalization group 
transformations explicitly is based on the averaging over rapidly oscillating degrees of free- 
dom. Background field method in the path integral framework || is usually associated with 
this procedure. The Hamiltonian counterpart of the background field method was also em- 
ployed in some problems, for example, in connection with gauge fields on a torus ]3], in 
soliton quantization || and also in variational calculations in gauge theories || . 

The averaging procedure is common in quantum mechanics with finite number of degrees 
of freedom. An instructive example |8| is the system with Hamiltonian 

H=\pl + \p 2 y + \x 2 y 2 - (1.1) 

The potential energy for this system is degenerate along the coordinate axes. Far from the 
origin, say, at \x\ ^> 1, the potential valley becomes very narrow and the wave function varies 
along the valley much slower than in the transverse direction. So, the variables are naturally 
separated in the slow and the fast ones. The last two terms in (|1 . 1|) is the Hamiltonian for the 
fast mode y: H h = \y 1 y + \x 1 y 1 . Its ground state energy induces an effective potential for the 
slow degree of freedom: V e s(x) = \x\/2. Thus, the degeneracy of the potential is lifted and 
the system with Hamiltonian in fact, has a discrete spectrum [0]. The above procedure 
can be generalized to a field theory ||, where the problem of lifting of the classical vacuum 
degeneracy also can be addressed || |§. In this paper we use the averaging procedure in a 
systematic way to renormalize Schrodinger equation in ip A and in the Yang-Mills theories. We 
consider pure gauge theory mostly for the sake of simplicity and the introduction of quarks 
should not cause any difficulties, since suitable methods to treat fermions in the functional 



Schrodinger picture are known [10 



It is worth mentioning that the averaging over fast modes in not the only way to renor- 
malize Schrodinger equation. As usual, the renormalization group transformations consist in 
a proper modification of the Hamiltonian for low-energy degrees of freedom after elimination 



of the fast modes So, the renormalized Hamiltonian acts in the smaller Hilbert space 
than the bare one and the renormalization can be considered as a projection on this smaller 
space, or a reduction of the Hamiltonian to the block-diagonal form. General methods for 
partial diagonalization [[12], [13|, as well as for projection [14| of Hamiltonians exist. These 
methods were tested on simple quantum-mechanical systems ]T3|, [H| and applied to some 
many-body |13|, |18 | and field-theoretical light-cone 0] and equal-time [17]| Hamiltonians. 
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2 Scalar field 



The Hamiltonian of ip 4 theory is 
H = J d 3 x 

where II is canonically conjugate to the field variable: 



2 2 y 1 2 24 . 



(2.1) 



b(^),n(j/)] = iS(x-y). 



(2.2) 

Schrodinger equation for the Hamiltonian ( |2.1| ) requires regularization, thus we assume that 
modes with momenta larger than A are somehow excluded (the concrete prescription will 
be given below). In order to exclude also the modes with momenta larger than /x, where 
A 3> \l ^> m, we can explicitly solve Schrodinger equation for the fields containing only these 
modes and then average the Hamiltonian with the wave function for the high-energy degrees 
of freedom. Denote by (p, H and by 0, 7r the slow and the fast variables - the field components 
which contain modes with momenta p < [i and /i < p < A, respectively. According to the 
conventional assumption of the background field formalism || , slowly varying fields (p satisfy 
classical equations of motion. 

Extracting the Hamiltonian for the high-energy degrees of freedom and expanding it to 
the quadratic order in <p (which is equivalent to the lowest order perturbation theory in A), 
we get: 



d x 



Since this Hamiltonian is quadratic, the ground state wave function has the form 



exp 



d 3 xd 3 y <f>(x)K(x, y)<f>(y) 



(2-3) 



(2.4) 



Substituting this expression in the Schrodinger equation 

d 3 x 



E h V h , 



(2-5) 



we obtain for the operator K: 



K 



1 \ 1/2 



the ground state energy being equal to 



E h 



Tr K. 



(2.6) 



(2.7) 



As in the quantum-mechanical example in the introduction, the ground state energy of 
large-momentum modes renormalizes the low-energy Hamiltonian. When \x is much larger 
than any other scale in the problem, only the ultraviolet divergent contributions are essential. 
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UV divergences can be easily extracted by the heat kernel method, which is based on the 
following representation for the square root of the operator: 

1 1 f°° dr _ tK 2_s1_ 



K = lim — / — e" TA "~ . 2.8 



\2e 27TV2 J r 3/2 

When we calculate the trace, the field-independent term merely renormalize the zero-point 
energy, 

r d 3 p y/p 2 + m 2 

Eo = jj^r — 2 — • 

Below this field-independent contribution is ignored. The remaining part of the effective 
potential can be calculated, for example, expanding 



V eS {(p) = E h = const - J -^ Tr exp -r (-d 2 + m 2 + ~ Ay? 2 ) 



(2.9) 



in A. This procedure involves evaluation of the momentum integrals, where the momentum 
cutoff can be introduced as a lower bound of integration over r. Therefore, the prescription to 
integrate over r from A~ 2 to fi~ 2 naturally accounts for the presence of only large-momentum 
modes in cf). 

In fact, direct expansion in A is not the most convenient way to extract the essential 
terms in the effective potential. The UV divergences are governed by the behavior of the 
integrand in fl2.9|) at small r and can be easily captured by DeWitt-Seeley expansion 0. In 
particular [ |19|| : 



(x\e-^- d2+v ^\x) 1 



l + Vr+ [\v 2 - ^d 2 V)r 2 + 0(T :i 



(4ttt) 3 / 2 

Substituting \(p 2 /2 for V, we get from eqs. ( |2.9|) , ([2.10|) : 



(2.10) 



1 r dr e ~ m2r r / 1 1 \ 

V eS {(p) = const -—^l-^ j—^ j Sx (- A^ 2 r + -A 2 ^r 2 j + 0(1//, 2 ), (2.11) 

The divergent terms lead to quadratic renormalization of the mass and logarithmic renor- 
malization of the coupling: 

Q\2 A2 

A - = A -32^V (2 ' 12) 
This result, of course, gives the correct expression for the j3-i unction in the ip 4 theory. 



3 Yang-Mills theory 

The Hamiltonian formulation of the Yang-Mills theory is most simple in the temporal gauge 
A Q = 0. Then the canonical variables are gauge potentials Af (x) and electric fields Ef(x) 
(we consider SU (N) gauge group, so B = 1, ... , iV 2 — 1): 

[A?(x),E?(y)] = i5 AB 8 ij 8(x-y). (3.1) 
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The Hamiltonian is 



H = J#z\^EtE? + ±F$F$), (3.2) 

where = d%A A — djA A + f ABC AfA^. The wave functions of physical states are also 
subject to the Gauss' law constraint: 

DiEf * = 0. (3.3) 

The covariant derivative Di acts in the adjoint representation: D AB = 5 AB di + f ACB Af . 

In the "coordinate" representation, the wave function is a functional of the gauge poten- 
tials and the electric field operators act as the variational derivatives: E A (x) = —i8/5Af(x). 
In order to find the dependence of the vacuum wave functional on the field modes with high 
momenta, we proceed in the same way as in the previous section. The slowly varying fields, 
Ai, are assumed to contain only modes with momentum p < p, and to satisfy classical 
equations of motion: 

DiF i:j = 0, (3.4) 

where Dj is the covariant derivative with respect to Ai and Fij is the corresponding field 
strength. It is convenient to rescale the fast variables cij, which contain the modes with 
momenta fi < p < A, so that Ai = Ai + gcii. To preserve the canonical commutation 
relations for the fast modes, the high-energy components of the electric fields should be 
rescaled by 1/g: Ei = Ei + - e«. 

We assume that the scale fi is sufficiently large, so that the running coupling g(/i) is 
small, and we can use perturbation theory for high-energy Hamiltonian. Expanding the 
Hamiltonian for rapidly oscillating variables in g to the leading order, we get: 



H h = / d'x 



\ 44 + \ 4(~D% - 2 Fij + DiDj) AB af 

Here F AB = [Dj, Dj] AB = j ACB Ff- acts as a matrix in the adjoint representation. 
Denote by L the quadratic form of the potential in eq. (|3.5|): 



(3.5) 



L{j = -D% - 21',, + DiD 3 , (3.6) 
Then the ground state wave functional for the fast variables is 

^^exp^-^aL 1 ^, (3.7) 

where summation over color and spatial indices and integration over spatial coordinates is 
implied. The wave functional ( |3.7|) satisfies Gauss' law constraint up to the two first orders 
in g: 

^ Di4 + DiEf + f ABC A B ef ) = 0. (3.8) 
The 0(1/(7) part of the Gauss' law generates transformations 



^ — > cii + DiUj. (3.9) 
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The invariance of the wave functional under these transformations follows from the identity 
LijDjU = ||. Terms of order g° in the Gauss' law generate usual gauge transformations 
of the background fields and homogeneous transformations of the fast variables: 

Ai -> fl\di + Ai)Q, ai ^VL ] aiVL. (3.10) 

The invariance of ( |3.7|) under these transformations is evident, since the kernel of the operator 
L 1 / 2 transforms homogeneously: L 1 / 2 (x, y) — ■> f2^(x)L 1//2 (x, y)Vt{y). To recover the invariance 
of the wave functional under the transformations generated by higher order terms in the 
Gauss' law, it is necessary to include higher orders in the Schrodinger equation. 

The effective potential generated by averaging of the Hamiltonian with the wave func- 
tional ([3.71) is 

Ves(A) = ^ Tr L 1/2 . (3.11) 

Although the operator L has a large number of gauge zero modes, they do not cause any 
problems with the calculation of the trace of its square root, hence, there is no reason to fix 
the gauge and to introduce ghosts as in the conventional background field method, where one 
usually deals with logarithms of the operators. The field- independent part of the effective 
potential corresponds to the zero-point energy and is not traced below. If the scale /i is 
sufficiently large, only UV divergent contributions are essential. The only contribution of 
this kind, the gauge coupling renormalization, can be easily extracted by the heat kernel 
method, as in Ref. |J, where a similar expression was considered in the context of the 
variational approach to the QCD vacuum. 

DeWitt-Seeley coefficients for the operator L were calculated in Ref. [2IJ, where a more 
general operator, 

1 



L -(a) = —D 5ij - 2F ij + (!--) A^, (3.12) 

was considered. Taking the limit a — > oo of the small r expansion of the heat kernel given 
in Ref. @]0, we find: 

— If (It f 

V eS (A) = const--— / — - d 3 x tr(x\e- rL \x) 



47r l/2 J r 3/2 

1 r dr 1 / o 11JV 



const + r 2 / JPxF£f£ + 0(r 3 

6 J 3 3 



4?r l/2 J r 3/2 ( 47rr )3/2 
1 1 ]\J A2 i r 

C ° nSt - W ln ^I / ^W + OCV/* 2 )- (3-13) 



Here tr denotes the trace with respect to color and to spatial indices. The averaging over 
the fast modes reproduces the usual coupling constant renormalization: 

1 1 11JV , A . _ N 

— = ^ln-, 3.14 

g 2 eS g 2 24vr 2 n 

as expected. 



*In fact, the first, field-independent, DeWitt-Seeley coefficient of the operator L(a) diverges as a — > oo. 
Although we do not trace the field-independent contribution, the origin of this divergence is clarified in 
Appendix. The correct prescription consists in simply dropping the divergent term. 
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4 Discussion 



The most common application of the background field techniques, that to QCD sum rules 
[ETJ, is based on the parametrization of the slowly varying fields by vacuum condensates 
[22|j . In the Hamiltonian picture, this would correspond to treating rapidly oscillating fields 



perturbatively, as above, and parameterizing the wave functional for the low-energy degrees 
of freedom by a finite number of parameters. Probably, this approach may provide some 
useful information. 

One may try to parametrize the vacuum in Yang-Mills theory by a finite number of 
parameters using some reasonable ansatz for the ground state wave function and then to 
find an approximate solution of the Schrodinger equation from the variational principle (for 
a recent discussion of the variational approach to the QCD vacuum, see p3| ||). Apparently, 
an accuracy of such kind of approximations is not under control. In particular, the correct 
renormalization of physical quantities is not automatically guaranteed in the variational 
calculations. Our results show that any approximate wave functional compatible with the 
asymptotic freedom must depend on the large-momentum modes of the fields in a very 
specific way - as given by eq. (|3[7|). This property substantially restricts valid functional 
form of variational ansatze for the ground state in Yang-Mills theory. 
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Appendix A Heat kernel of the operator L 

In order to relate the operators L and L(a), it is convenient to introduce the non-Abelian 



transverse projector 24 



Sy-DirLDj. (A.l) 



G ij = -D 2 5 ij -2F ij . (A. 2) 



Denote L(l) by G: 

r 

Then, with the use of the equations of motion ( |3.4| ) and the commutation relations [Dj, Dj] = 
Fij, the following equality can be proved: 

Lij = GikPkj = PikGkj, (A. 3) 

So, the operator L is a transverse projection of the operator G. 

As a consequence of eq. ( |A.3| ), and since P is the projection operator, 



e - tL = e ~ tG P. (A.4) 
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One the other hand, eqs. 
consequently, 



( pOD and flOp imply that A^j 



Lij(a) — GikPkj + — Gik(5kj — Pkj)- 



kj 



Since P and 1 — P are the orthogonal projectors, 



-tL(q) 



e ~ tG P + e 



- — G i 



P kj ) and, 



(A.5) 



(A.6) 



The diagonal matrix elements of the second term on the right hand side of eq. (|A.6| ) can 
be expanded in local operators. By dimensional reasons the parameter of this expansion is 
r/a: 



tr (x\e-« G (l -P)\x) = 



n=0 



OJx) 



(A.7) 



where tr denotes the trace with respect to spatial and to color indices and O n (x) are gauge- 
invariant operators of dimension 2n. As there are no such operators of dimension 2, only the 
first term proportional to a 3 / 2 survives the limit a — > oo. It is field- independent and can be 
calculated in the momentum representation. Finally, we get: 



tr(x|e tL \x) 



lim 

ct— >oo 



tr (x\ e 



-tL(q) 



x) 



(N 2 - l)a 3 / 2 
(4ttt) 3 / 2 



(A.J 



The small r expansion of the right hand side is given in Ref. |20[ up to the two first terms 
and, indeed, has a finite a — > oo limit. 
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